Luttinger Liquid in the Core of Screw Dislocation in Helium-4 
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On the basis of first-principle Monte Carlo simulations we find that the screw dislocation along 
the hexagonal axis of an hep 4 He crystal features a superfluid (at T — > 0) core. This is the first 
example of a regular quasi-one-dimensional supersolid - the phase featuring both translational and 
superfluid orders, and one of the cleanest cases of a Luttinger-liquid system. In contrast, the same 
type of screw dislocation in solid H2 is insulating. 

PACS numbers: 75.10.Jm, 05.30.Jp, 67.40.Kh, 74.25.Dw 



(N 



5h 



1 

C 

o 
o 



(N 
> 

VO 
0\ 
(N 

in 
O 

o 



X 



The remarkable observations of a non-classical moment 
of inertia in solid 4 He [ESS III! have sparked a wave 
of interest in the supersolid state of matter (see review 
@ and references therein). According to the standard 
paradigm Q , the uniform supersolid state develops if the 
energy gap for creating vacancies vanishes in the ground 
state of a crystal, leading to a finite concentration of sta- 
ble and delocalized vacancies. In the past, however, ex- 
periments did not find such zero-point vacancies toL Re- 
cent Monte Carlo simulations of solid 4 He [1, Holllll. HH 
have further limited possible mechanisms for a uniform 
supersolid state: A perfect hep 4 He crystal is an insula- 
tor with a large energy gap toward vacancy creation and 
an even larger one for interstitials. When vacancies are 
introduced by hand into a perfect hep crystal, they un- 
dergo phase separation [9( . This provides strong evidence 
that a 4 He crystal does not conform to any standard su- 
persolid scenario. 

The situation is quite different for defect-ridden, or 
polycrystalline hep 4 He. Simulations showed that grain 
boundaries are generically superfluid [l3| , while a recent 
experiment observed superflow in non-uniform samples 
and attributed it to the presence of grain boundaries [lj] . 
In general, two types of supersolid may exist - regular 
and glassy. While in both cases translational symmetry 
is broken, regular supersolid is characterized by two or- 
ders - discrete translational and superfluid. Thus, due to 
lack of discrete translational order, the superfluid grain 
boundaries do not constitute an example of a regular su- 
persolid. They should rather be considered as superglass. 

One may conjecture that (straight) dislocations, fea- 
turing translational order, might host SF along their 
cores [TH, [l|| . More recently, a phenomenological model 
was proposed, where crystal deformations pV7| . particu- 
larly inside a dislocation core promote SF see also 
Ref. [19( on role of 3 He in inducing superfluidity in dis- 
location cores. None of the previous theories, however, 
was quantitative and conclusive. The model [l?], 13 con- 
trasts with our observation of insulating grain boundaries 
with small tilt angles because this type of boundary can 
be considered as a wall of rarely spaced edge dislocations 



[131 [20J] . We are not aware of any suggestion that screw 
dislocations have any advantage over edge dilocations in 
terms of their SF properties; on the contrary, they were 
often excluded from the list. Thus, the question remains: 
Do extended supersolid defects characterised by transla- 
tional and superfluid orders exist in 4 He crystals? In this 
Letter, we give an affirmative answer to this question, 
i.e., we provide theoretical evidence of SF inside the core 
of a screw dislocation aligned with the c— symmetry axis 
of hep 4 He. 

Our grand-canonical Monte Carlo simulations are 
based on the worm algorithm, i.e., a Path Integral Monte 
Carlo technique defined on the configuration space for the 
single-particle Matsubara Green function [2l| . For spa- 
tial imaging, we employ the winding-cycle maps, i.e., av- 
eraged set of points measuring instant spatial position of 
a randomly selected element (bead) of the worldline with 
non-zero winding number (winding cycle). The wind- 
ing cycles define the superfluid response [22j], so that 
the winding-cycle map allows one to visualize the spa- 
tial distribution of the superfluid component. Qualita- 
tively speaking, a winding-cycle map represents the local 
superfluid response. 

Sample geometry. The geometry of our sample is il- 
lustrated in Figs. Q]and[2] The sample is periodic in the 
z-direction, with 8 pairs of basal planes. In order to rule 
out finite-size effects, we have simulated samples with 3 
and 4 pairs of basal planes as well. 

Cell periodicity in the xy plane is incompatible with 
the presence of a single dislocation. One could restore 
it by simulating a dislocation pair, but we opt instead 
for a single dislocation, using a non-periodic cell in the 
xy plane, as this allows us to a) simulate samples with 
the rotational symmetry of a hep crystal b) obtain robust 
results with smaller samples. 

Our simulation sample consists of an ideal hep crys- 
tal with a screw dislocation in its center along the c 
axis. Specifically, we, first, superimpose two identical 
triangular layers separated by half of the c-axis period 
a z = ^/8/3a along the c-axis and shifted by (a/2, a/y/3) 
in the xy-plane, a being the lattice period along the x- 
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FIG. 1: (Color online. For the best perception, look from 
a distance.) Columnar winding-cycle map (blue dots) in the 
core of the screw dislocation in solid 4 He at T — 0.25 K and 
density 0.0287 A -3 . View is along the hexagonal axis, with 
the core parallel to it. Shown with red dots (in the lower 
half of the plot only) are the atomic positions in the initial 
configuration. At distances ~ 3a from the core, the atomic 
positions averaged over the imaginary time are only slightly 
shifted against those of the original configuration. The unit of 
length is 1 A. 



axis. The crystal, then, has been cut along the rrz-half- 
plane, whose edge (located at the center of one of the 
hexagons formed when all basal layers are projected to 
a single xy-plane) is the dislocation core. Atoms have 
been displaced so that, upon completing a full revolu- 
tion around the core, each atom advances by a lattice 
period a z . This procedure creates an ideal (classical) 
screw dislocation in essentially infinite medium. Then, 
all particles located outside a (pencil shape) cylinder, 
with the dislocation core being its axis of symmetry, 
(and inside the rectangular simulation cell) have been 
pinned. These particles are not moved in our Monte 
Carlo simulation, but they do interact with the rest of 
the system via the helium pair potential. The purpose of 
these "frozen" particles is that of confining our sample 
to the inner (cylindrical) part of the cell, while, practi- 
cally, completely eliminating the effects of open boundary 
conditions. Our rectangular simulation cell is twice the 
size of the physical sample inside the cylinder. The ini- 
tial numbers of (physical) atoms used in simulations are 
N ini = 384, 768, 1912. 

Superfluid response. The simulations were performed 
at the number density 0.0287 A -3 , that is, in the close 
vicinity of the melting point. The lattice period along 
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FIG. 2: (Color online.) Columnar winding-circle map in the 
core of screw dislocation. View is along the x-axis in the 
basal plane — perpendicular to the core. Shown with red dots 
(in the lower half of the plot only) are the atomic positions in 
initial configuration. The unit of length is 1 A. 



the ^-direction is then a — 3.666 A. The convenience of 
working at the melting point is that we know the chemical 
potential, /i = 0.02 K, from the simulations of the liquid 
phase at its freezing density [2lj . Our simulations span 
a temperature range from 0.2 to 1 K. 

In Figs. [1] and [5] we present the projected winding- 
cycle maps on the xy- and yz-planes, together with the 
initial configuration of the atoms. We see that the super- 
fluid density in the dislocation core is strongly correlated 
with the insulating environment, respecting the hexag- 
onal symmetry modified by the presence of the screw 
dislocation. The superfluid density is most robust along 
the bonds of the ideal structure. 

Quantitatively, our data can be analyzed using the con- 
cepts of Luttinger liquid theory. Luttinger liquids are 
characterized by two parameters, the superfluid stiffness 
A s = n{ 1D ^ /to (where is the ID superfluid den- 

sity and to is the helium atom mass) and the compress- 
ibility k. Both can be extracted from the distributions 
Pw(W), Pn(N) of the winding W (along the core) and 
particle number N fluctuations directly obtained in the 
simulation: 



Pw(W) ~ exp 
Pat (TV) ~ exp 



(3N 2 
2Lk 
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where N is taken relative to its mean (N) [23], L is 
the sample length along the c-axis and /3 stands for in- 
verse temperature. We fitted the collected histograms 
Pw(W), Pjy(N) with Gaussians and have obtained the 
superfluid stiffness and the compressibility. The Lut- 
tinger parameter is given as = tt^/Ak. Our data (for 
the different temperatures, and different system sizes in 
the z-direction) yield Kl = 0.205(20), where the finite- 
size (along the c-axis) and statistical errors are combined. 
The initial number of (updatable) particles in the state 
shown in Figs. [T] and [5] was N — 768 corresponding to a 
relatively tight confinement in the icy-plane. We clearly 
see superfluidity in samples with much larger number of 
particles in the xy-planc (1 91 2 updatable particles in the 
sample with 3 pairs of basal planes). In particular, the 
superfluid stiffness was found to be the same within the 
error bars (of about 20%). However, the accuracy of data 
for k was not sufficient for resolving finite-size effects in 
Kl due to the icy-confinement. 

For the value of the chemical potential corresponding 
to the melting point, the average equilibrium number of 
particles in the 8 basal-pair sample was found to be about 
770. By varying chemical potential and, in particular, 
lowering it below the point where N = 768, we saw no 
qualitative difference in the superfluid response, which 
is fully consistent with the fact that Kl < 1/2 (in the 
Luttinger liquid, an external periodic potential is rele- 
vant for opening the insulating gap only at Kl > 1/2 
[H, [H, Ha [13). A sample half the size in z-direction 
yielded the same qualitative conclusion. In the 3 basal- 
pair sample with much larger number of particles in the 
icy-plane, we observe that the core dopes itself with va- 
cancies (of the order of 3-4) , rather than with interstitial 
atoms, apparently due to lower effective chemical poten- 
tial in the core. We do not attribute special significance 
to this fact because by increasing the chemical potential 
by 4 K, far smaller than the interstitial energy in the 
hep solid, the system can be made to have the original 
(classical solid) particle number without changing the su- 
perfluid nature of the state. 

From our data for A s we deduce that ni 1D ^ is about 
1 particle per 1 A, which is equivalent to saying that su- 
perfluid phase in the dislocation core involves nearly all 
atoms in a "tube" of diameter 6 A. More quantitatively, 
one can assign the superfluid response to partial con- 
tributions from the coordination shells. Nearly 90% of 
points in the winding-cycle map in Fig. [T] belong to the 
first three coordination shells with the following division: 
35 % in the first shell, 29 % in the second one, and 26 % 
in the third one. 

Special care should be taken to make sure that the 
observed superfluid response is not an artifact of poor 
thermalization of the worldline configurations in the dis- 
location core. This concern stems from our previous ex- 
perience with edge dislocations, where we saw that a 
superfluid glassy region in the core — created either by 
hand, or as a result of quenched relaxation of the initial 
configuration — remains essentially metastable, at least at 



temperatures ~ 0.2 K. To rule out such a scenario, we 
first relaxed the initial configuration by treating atoms as 
classical particles with a repulsive potential ~ 1/r 6 , and 
minimizing the energy of the configuration. [In the case 
of edge dislocations, this protocol leads to an insulating 
groundstate.] Our observation was that in both cases — 
with or without classical relaxation — the Monte Carlo 
process converges to one and the same result. Another 
indication that our superfluid signal is not an artifact 
of poor thermalization comes from the perfect regularity 
of the superfluid density map, in combination with its 
strong correlation to the insulating environment. This 
outcome rules out concerns based on a lack of thermal- 
ization in the glassy state. Analogous simulations have 
been performed for the hep para-hydrogen crystal, using 
the same sample preparation procedure, at a density of 
0.0261 A~ 3 , corresponding to the T — » equilibrium den- 
sity. In this case, the dislocation core is found to be an 
insulator doped with vacancies (forming a commensurate 
density wave). 

Discussion. The observation of a superfluid dislocation 
in a 4 He crystal is an interesting example of a regular 
quasi-one-dimcnsional supersolid. However, the question 
is whether this observation is relevant to the effect of 
non-classical inertia [l[ . Superfluid dislocations can be a 
part of the disordered superfluid network, together with 
superfluid grain boundaries, ridges, and liquid (or glassy, 
at higher pressure) pockets, etc. Considering a network 
of dislocations only, the effective three-dimensional su- 
perfluid fraction of such a network is limited by the vol- 
ume fraction of the atoms located at the cores of dislo- 
cations. This fraction is likely to be significantly smaller 
than ~ 1 %. The dislocation density in solid 4 He typ- 
ically varies over a range as wide as 10 6 10 10 cm -2 
H, [2^, [13, [Ml, [H[ . It is, however, also pos sible to grow 
crystal without screw dislocations at all |33l | . Taking the 



ID superfluid density of n. 



(ID) 



1 A 1 observed in our 



simulations, we estimate the required order-of-magnitudc 



dislocation density ~ l/l 



3 x 10 12 cm 2 in order 



to account for nf /n = (n[ 1D> nd)/n ~ 1% of the su- 
perfluid fraction observed in [l|. Here I is the typical 
distance between dislocations. 

The Luttinger liquid parameter is all we need to char- 
acterize superfluid properties of the dislocation network 
consisting of segments of length I ^S> a. As it has been 
pointed out by Shevchenko 3J] , in this system one has to 
distinguish between static and dynamic responses. Below 
the thermodynamic transition temperature T c ~ T*a/l 
(which can be also written as T c ~ T* [ni 3D ' ) /n] 1 / 2 , in 
terms of the zero-point superfluid fraction ris D ^/n) the 
network features a non-zero static superfluid response 
[34] ]. Here, T* is the characteristic temperature of the 
superfluid helium liquid. In our simulations we find that 
the core has robust phase coherence properties at all tem- 
peratures below 1 K, which allows us to set T, ~ IK (the 
continuation of the A-line in liquid helium to the region 
of solid densities gives T\ w 1.5 K [13|). This estimate is 
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in agreement with the small value of Kl . 

In the temperature range T c < T < which is very 
broad for large aspect ratio I /a, the state of the sys- 
tem can be viewed as a vortex tangle with 'vortex cores' 
(which are not well defined geometric objects because 
of the natural uncertainty ~ I in their sizes and posi- 
tions) pinned in the solid bulk. We refer to this state as 
Shevchenko state (34J. In contrast to the original purely 
classical analysis of Ref. [34| , we emphasize that the dy- 
namics of Shevchenko state should be driven by discrete 
quantum-tunneling phase slip events (instantons), with 
the typical relaxation time r ~ (T./T) 8 /^ -1 /T*, as it 
follows from the microscopic analysis [35| . In the low- 
frequency limit lot <C 1 the state is essentially normal. 
At "high" frequencies lot 3> 1 the response is nearly in- 
distinguishable from the standard superfluid [34| . Note 
also, that for Kl <C 1 and T<T, the phase-slip time can 
easily exceed any realistic experimental time-scale mak- 
ing it very hard to detect the genuine thermodynamic 
transition at T c . 

While being focused on one isolated screw dislocation, 
we note that small twist-angle grain boundary can be rep- 
resented as an array of screw dislocations 20| and, there- 



fore, must also be superfluid with an additional feature 
of tunneling coupling between the cores. Such bound- 
aries might be the most realistic objects for studies of 
superfluidity of screw dislocations. 

In summary, we have found that the screw dislocation 
in hep solid 4 He supports (at T — > 0) superfluid transport 
of 4 He atoms along its core, yielding thus the first realistic 
example of a regular continuous supersolid. Analogous 
dislocation in the hep solid H2 is found to be an insulator. 
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